Strong anisotropy in surface kinetic roughening: analysis and experiments 
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We report an experimental assessment of surface kinetic roughening properties that are anisotropic 
in space. Working for two specific instances of silicon surfaces irradiated by ion-beam sputtering 
under diverse conditions (with and without concurrent metallic impurity codeposition), we verify 
the predictions and consistency of a recently proposed scaling Ansatz for surface observables like the 
two-dimensional (2D) height Power Spectral Density (PSD). In contrast with other formulations, this 
Ansatz is naturally tailored to the study of two-dimensional surfaces, and allows to readily explore 
the implications of anisotropic scaling for other observables, such as real-space correlation functions 
and PSD functions for ID profiles of the surface. Our results confirm that there are indeed actual 
experimental systems whose kinetic roughening is strongly anisotropic, as consistently described by 
this scaling analysis. In the light of our work, some types of experimental measurements are seen 
to be more affected by issues like finite space resolution effects, etc. that may hinder a clear-cut 
assessment of strongly anisotropic scaling in the present and other practical contexts. 
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Surface kinetic roughening, namely, the lack of charac- 
teristic length and time scales in the fluctuations of the 
values of an interface topography, is quite a pervasive 
phenomenon, from hard to soft-condensed matter physics 
systems. It occurs for instance in thin-film production^ 
crack formation,^ fluid dynamics,'^ or growth of cellu- 
lar aggregates.'* Naturally, in each specific context, the 
meaning and even the dimensionality of the relevant sur- 
face or interface needs to be unambiguously defined. In 
our case, we will consider as reference systems proper 
two-dimensional surfaces of e.g. thin films, produced by 
chemical or vapor deposition techniques, epitaxy, ion- 
beam erosion, etci^ 

In spite of the large efforts devoted to the understand- 
ing of kinetic roughening, a number of issues remain 
poorly understood. This lack of complete knowledge is 
perhaps one important reason for the relatively moder- 
ate experimental impact that this general phenomenon 
has met, see \^ and references therein. One such as- 
pect is the phenomenon of anisotropic scaling, namely, 
when the properties associated with scale invariance in 
the fluctuations of the surface height change with the sub- 
strate direction one is looking along. Technically, this is 
reflected in the fact that the ensuing power-law behav- 
ior of observables in the system^ are characterized by 
direction-dependent critical exponents. 

Perhaps the most obvious class of systems in which 
topographical anisotropics occur correspond to those in 
which morphological instabilities induce different typi- 
cal length-scales along different directions. For instance, 
for sofid targets eroded by ion-beam sputtering (IBS), 
an oblique ion flux induces ripple formation^i^ breaking 
the symmetry between the substrate directions along the 



ripples and perpendicular to them. Note, however, that 
for such pattern-forming systems scale-invariance does 
not occur at length scales that are comparable with the 
typical wavelength of the pattern. Nevertheless, kinetic 
roughening can still occur simultaneously at length scales 
that are, rather, well separated from such a preferred 
scale. Examples are known in e.g. the context of the 
Kuramoto-Sivashinsky (KS) equation;^ frequently used 
to model IBS systems;^ under such conditions, it is just 
natural to expect the ensuing kinetic roughening proper- 
ties to take on an anisotropic formi^^ 

For morphologically stable systems, one can con- 
ceive of anisotropic scaling behavior due to e.g. physical 
anisotropics in the external flux that is driving the sur- 
face dynamics. This is the case, for instance, for macro- 
scopic geological systems, in which an external matter 
flux has a direction that is preferred e.g. by gravity on 
a tilted landscape?ii Anisotropics can also exist in the 
relaxation mechanisms that are intrinsic to the system: 
e.g. anisotropics in surface tension in phenomena of crys- 
tal growth from a meltji^ or in barriers to relaxation by 
surface diffusion, even for a negligible (unstable) Ehrlich- 
Schwoebel effect, in the case of epitaxial growthii 

Given the above, remarkably few reports exist in which 
anisotropic kinetic roughening has been unambiguously 
assessed (see e.g. [ij] and references therein), specially as 
taking place in the steady state of the system, which is 
what we will call strong anisotropy (SA). From the the- 
oretical point of view, this is the most interesting man- 
ifestation of anisotropic behavior, since it corresponds 
to the behavior of a system in the thermodynamic (or, 
for out-of-equilibrium systems, more properly "hydrody- 
namic") limit. Perhaps one of the complications in the 
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identification of SA is the fact that this type of behav- 
ior had been previously encoded in a scahng Ansatai^ 
that originates in the study of critical dynamics of equi- 
librium statistical-mechanical systemsii^ Such a formu- 
lation is quite powerful from the theoretical point of view 
(enabling analysis of scaling properties for arbitrary sub- 
strate dimension, etc.) but is not particularly natural for 
the characterization of actual two-dimensional surfaces. 
In this way, the experimental thin-film, or surface-science 
community does not have a clear guiding principle that 
allows to identify correctly the behavior that should be 
expected for each observable under conditions for strong 
anisotropy. Comparison with theoretical models is more- 
over hampered. A clear-cut picture is particularly re- 
quired in situations in which finite size or analogous ef- 
fects may interfere with an underlying scaling behavior. 

Some formulations of SA adapted better to the thin 
film context are available,^ but have been tested onto 
simplified models only,-^ so that their generality has not 
been checked, nor has a systematic study been done on 
the relation between different observables within such a 
framework. Specifically, given the widespread current use 
of scanning probe microscopies (AFM, STM) and of X- 
ray diffraction techniques, it is important to understand 
the relationship between the scaling behaviors of the two- 
dimensional Power Spectral Density (PSD) of the surface 
height and the PSD of one-dimensional cuts along dif- 
ferent substrate directions, and their relation with real- 
space height correlation functions (see below for precise 
definitions). This is due to the fact that many experi- 
mental setups are often designed to measure either the 
former (X-rays), or the latter (AFM/STM)^ Therefore, 
a fully consistent analysis and a quantitative compari- 
son with theoretical models cannot be done unless this 
relationship is fully clarified. 

The purpose of this work is two-fold: first, consider- 
ing a recently proposed anisotropic scaling Ansatz that 
has been tested against linear and nonlinear models, 
we elucidate the behavior that is expected for one- and 
two-dimensional PSDs in a shape that is natural to thin- 
film analysis techniques, and the relation with the be- 
havior of real-space height correlation functions. Note 
that anisotropic scaling is a physical property that exper- 
imental systems may or may not have, independently in 
principle of comparison with theoretical models (which, 
in turn, may or may not have it on their own). 

Once this framework for analysis is clarified, we ad- 
dress the experimental observation of SA in two cases. 
They both correspond to experiments of surface erosion 
of silicon targets by ion-beam sputtering. But one of 
them one is within a morphologically stable condition 
in which metallic contaminants are suppressed, while a 
second one corresponds to erosion with concurrent de- 
position of metallic contaminants leading to a morpho- 
logical instability. At any rate, the occurrence or not of 
the instability will not play any role in our discussion, 
as we will focus on length scales at which scale invari- 
ance holds in each case. Thus, we are able to discuss in 



detail various aspects on the occurrence of SA that may 
have hampered a more clear assessment of this property 
in the past. Based on these results, we propose forms of 
analysis that seem less susceptible for complications and 
ambiguities. 

This paper is organized as follows. In Section HI we 
discuss the anisotropic scaling Ansatz that will be em- 
ployed for data analysis in later parts of the work. We 
describe its implications for height measurements that 
are performed both in real and in Fourier space, for the 
full two-dimensional surface, and/or for line profiles of it. 
Section [n] is then devoted to the application of this scal- 
ing Ansatz to the description of our experimental data. 
Finally, we present a discussion of these results from a 
general point of view and conclude in Section IIIIl Sup- 
plemental Material is provided^* that contains morpho- 
logical data collected by AFM and STM as well as some 
auxiliary plots, that has been employed in order to pro- 
duce the graphs and the discussion provided in Section 
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I. ANISOTROPIC SCALING ANSATZ 

Consider a surface whose height above a reference 
plane is given by a height function h{r), where r — (x, y) 
is a point on a reference substrate plane. If the surface 
morphology is disordered;^ a quantity that conveniently 
characterizes the fluctuations of the height around its av- 
erage value is the height-difference correlation function 

G(r) = ([Mr + ro)-Mro)]'), (1) 

which is usually a function of relative distance, r = |r|, 
only. Here, brackets denote average over, say, different 
experiments carried out under the exact same conditions, 
and To — (a;o,j/o) is an arbitrary position on the sub- 
strate. Isotropic kinetic roughening corresponds to a sit- 
uation in which this correlation function grows as a power 
law with distance, G(r) ^ r^", where a is the so-called 
roughness exponent 4. This exponent has a direct relation 
with e.g. the fractal dimension of the surface, thus char- 
acterizing its asperity or roughness properties. Note, the 
term "kinetic roughening" suggests in particular that our 
surface is evolving in time. Nevertheless, unless stated 
otherwise, we assume in this Section that a steady state 
has been reached in which morphological properties do 
not change on average. 

For an anisotropic surface, the value of the roughness 
exponent changes with the direction along the substrate 
plane, thus one has, Gx{x) ~ x"^"" and Gy{y) ~ y^"", 
where 

Gx{x) = {[h{xQ + X, yo) - h{xo, yo)f), (2) 

and Gy{y) is defined analogously. The system is said 
to display Strong Anisotropy (SA) if indeed ^ ay, 
whereas we talk about Weafc Anisotropy (WA), when the 
steady state of the system is actually isotropic. 



3 



In our analysis we focus on the power spectral density 
(PSD) of the surface height, 



5(k) = (iMk)r), 



(3) 



where /i(k) is the space Fourier transform of h{r) — h, 
with h being the space average of the height, and k = 
(kx, ky) is wave- vector. Our aim is to formulate a scaling 
form for S'(k) that is "typical" of strongly anisotropic 
systems, and study how such a behavior reflects in other 
observable functions related with the PSD and/or the 
height-difference correlation function. For an strongly 
anisotropic system, we postulat&ii that the stationary 
PSD scales with wave-vector components k^ and ky as 



S [kx , ky 



2q. 



(4) 



where we refer to dix and Oiy as roughness exponents 
in momentum space, and u is a mere constant. Note 
that the "asymptotic to" (^) sign expresses an "equality" 
that holds, up to numerical constants, in an appropriate 
small-A; = |k| (equivalently, large-r) approximation, as 
customary in critical behaviorJ^ Note also that, in prin- 
ciple, the behavior expressed by Eq. (g]) can be directly 
checked from, say, AFM or STM data treated in Fourier 
space. Moreover, the PSD function S'(k) does behave as 
described by Eq. ([4|) for a number of representative lin- 
ear and non-linear models of anisotropic surfaces;^ As an 
example, consider a surface that relaxes by surface ten- 
sion along the x direction and by surface-diffusion along 
the y direction. The simplest evolution equation for the 
surface height that incorporates fluctuations in, say, an 
external flux of adatoms read a^^'^^ 



dth = Vxdlh 



(5) 



where h'x,y are positive constants related with surface ten- 
sion and surface diffusivity,^'^ and 77 is a noise term. For 
this linear model, it can be readily shown that, in the 
long time f — 00 limit, the PSD function takes exactly 
the form given by Eq. ^ with a^; = 1 and fiy = 2; more 
elaborate models also comply with the same functional 
form, with different values for the cix.y exponents, see 
[13 and Sec. HI] below. 

In many experimental systems, the two-dimensional 
PSD is too noisy to get a reliable estimation of the crit- 
ical exponents, and a different approach is preferred. 
A simple way to increase the signal-to-noise ratio is to 
take one-dimensional profiles of the surface topography 
h(x, y) along the x (or the y) direction and Fourier trans- 
form them. Given a point yo (or, equivalently, xq) one 
can show that the continuum limit of the PSD of this 
one-dimensional cut is obtained by integrating the two- 



dimensional PSD over the ky directioiiii 



Sx{kx) — / dky Sikxj ^y)- 



(6) 



Obviously, the same relation holds for the PSD of the 
one-dimensional cut along the y direction, i.e. Sy{ky). 



Inserting Eq. Q into ([6]), we obtain scaling laws for Sx 



and S*,, 



Sx(kx) 

Sy (ky ) 



-(2a^-C) 



'1/0 



(7) 

(8) 



In these relations, we have introduced the exponent ratio 
C = OLx/oiy that measures the degree of anisotropy of the 
surface, and will be referred to as an anisotropy exponent. 
We say that the system displays SA whenever C ^ 1, 
while C = 1 for isotropic systems. As shown belowfiS 
Eqs. d?]), (m are equivalent to 



Sx ij^x) ^ X + ) ^ 



Sy{kyl 



(9) 
(10) 



that provide the natural generalization to the SA case 
of the scaling behavior of the PSD of ID cuts of the 
surface in the isotropic case, in which ax = ay — a 

and Sx^y ^ kZjy"'^^^'^ Still, we have to prove that 
the anisotropy condition C 7^ 1 indeed corresponds to 
c^x 7^ Q^j/j for which we need to relate these real-space ex- 
ponents with their momentum-space counterparts, ax,y- 
The required relation is readily obtained from the defi- 
nition of the height-difference correlation function along 
one direction, e.g. along the x coordinate, Gx{x) above. 
To increase the precision of the estimate of Gx , this func- 
tion will be averaged over all values of the vertical co- 
ordinate, yo- The same procedure applied along the y 
direction gives Gy{y). From the exact relation between 
Gx and Sx^ 

2 f°° 

Gx{x) = - / dkx [1- cos(kxx)] Sx{kx), (11) 



we get the scaling law for Gx as a function of the rough- 
ness exponent defined in Fourier space 



Gx{x) X 



2Qx-C-1 



(12) 



Finally, by equating this last result with the scaling 
Ansatz for Gx in real space, we are able to relate the 
different roughness exponents 



2ax = "iux - C - 1, 



(13) 



and, by applying the same reasoning to the y direction, 
we get the relation 



2ay = 2ay - 1/C- 1. 



(14) 



Eqs. in]) and HI]) immediately convert ©, ^ into © 
and (jlOp. Moreover, it is simple to show that, in fact. 



^ ^x ^x 



(15) 



The anisotropic scaling Ansatz expressed by Eq. ^ 
can be shown to be equivalent to other that are known 
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in the literature on critical systems, see e.g. [ll|, [3], 
and references therein. Actually, it is an analogue for 
non-conserved noise fluctuations of the behavior found 
in so-called driven-difFusive-systenis»22, Nevertheless, an 
advantage of the present formulation is that one can re- 
late the scaling behaviors of different 2D and ID corre- 
lation functions with one another in a way that readily 
generalizes the isotropic behavior. In practice, this is 
a convenient consistency check that enables safe conclu- 
sions on the occurrence of the interesting property of SA 
in a given experiment, once compatibility ensues between 
the expected scaling behavior of the various observables 
as measured on the same set of data. We are not aware 
of such a type of analysis in the literature (see a related 
analysis of SA behavior in 1^ ) , so that the next section 
is devoted to the discussion of a couple of experimental 
systems in which SA scaling, as described in the present 
section, is indeed shown to occur. 

Note also that Eq. (|1]) has been postulated once the 
PSD function 5'(k) has reached a steady state and is 
time independent. The experimental data we report 
about in the next section satisfy this condition. Nev- 
ertheless, in order for this type of behavior to prop- 
erly become an anisotropic generalization of the so-called 
Family- Vicsek (FV) Ansatz that applies to isotropic ki- 
netic roughening,'* time dependence should be allowed 
for. The FV Ansatz is typically formulated in terms of 
the short and long time behavior for the surface rough- 
ness W^it) = {{h - hf) = J S{k) dk. Thus^ W tf^ 
for t < t^/^, while W ^ L" fov t :^ t^/^, where z is an 
independent exponent, i^/^ is an estimate of a length- 
scale below which non-trivial correlations have built up 
among height values at different substrate positions, and 
j3 — a/z. In the SA case, it is natural to consider cor- 
responding roughness functions Wx,y for one-dimensional 
cuts of the surface, so that one would expect Wx,y ^ t^"-^ 
for t < tVzx.H ^ while W^^y ~ i""" for t > iV^x,,; ^ xhis is 
indeed the behavior compatible with the steady state Eq. 
(j4|), for which one can show^^ that moreover Zy — Zx/C,, 
and thus fix = ctx/zx = Py = ay/ Zy. Hence, overall there 
are only three independent critical exponents character- 
izing a time-dependent SA surface, e.g. ax, C: and z^, 
from which all other exponents described in this section 
can be obtained. 



our case, these properties will correlate with the simulta- 
neous codeposition of metallic impurities, or the absence 
thereof, as recently shown.i^LS^ At any rate and as men- 
tioned above, even in the morphologically unstable case, 
it may be possible to find a range of scales within which 
scale-invariance (i.e., kinetic roughening) occurs, that is 
the one we will address. This is the case for large classes 
of pattern- forming systems.^ 

In each section below, the analysis of the experimental 
data, and the verification of the scaling Ansatz presented 
in Section HI has been carried out following the next sev- 
eral steps: 

• For each experimental condition the irradiation 
dose guarantees that the surface has reached a 
steady state (as reflected in the time independence 
of the correlation functions measured) . 

• A single sample is scanned (by STM or AFM) using 
different resolutions, and over regions of different 
sizes, which we will call windows. These scans pro- 
vide us with different images, encoded in matrices, 
of the same sample. 

• For each matrix, we compute two quantities: the 
two-dimensional PSD S{kx,ky), and the PSDs of 
one-dimensional cuts of the surface, Sx,y{kx,y)- The 
2D PSD is then averaged over all the images we 
have. For the ID PSDs, as explained in Section 
U we also perform averages over all the lines and 
columns of each matrix, respectively. 

• Next, we consider the projections S{kx,Q) and 
S{0,ky) of the two-dimensional PSD onto the kx 
and ky axes. According to our Ansatz, these pro- 
files should scale as power laws of kx^y with expo- 
nents 2ax and 2ay , respectively, that are estimated 
by fits over appropriate scaling ranges. 

• Using formulas (IT5|) and ([H]), we compute real- 
space roughness exponents ax,y and, through equa- 
tions © and (ITU|) . we predict the scaling behavior 
for the one-dimensional PSDs. 

• Finally, we verify if the one-dimensional PSDs in- 
deed scale in a form that is consistent with the two- 
dimensional PSD, as required by our scaling theory. 



II. EXPERIMENTS 

In this section we show that the previous scaling analy- 
sis does provide a consistent description of the anisotropic 
kinetic roughening properties of actual experimental sur- 
faces. In particular, we choose to produce the lat- 
ter through erosion of silicon targets by IBS at low- 
intermediate ion energiesi^ By tuning experimental con- 
ditions such as the angle of incidence of the ion beam 
onto the surface, average ion energy, etc. it is possible 
to obtain surfaces with varying topographical proper- 
ties, from disordered and rough to nanopatterned. In 



A. Morphologically stable case 

For the morphologically stable condition that we will 
study, the experiments were performed in a scanning tun- 
neling microscopy (STM) apparatus with a base pressure 
below 1 X 10~^° mbarr^^ The system is equipped with a 
differentially pumped fine focus ion source. The fine fo- 
cus ion beam exposed only the sample and thus impurity 
effects were avoided. A Si(lOO) sample was irradiated at 
room temperature with 2 keV Kr"*" ions at an ion inci- 
dence angle 9 = 81° with the surface normal. The repro- 
ducibility of the angle is better than 0.5°, the error of the 
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FIG. 1. STM top views of the surface morphology for the 
contamination- free sample. Image size is 1.3 x 1.3 nm^. 

angle's absolute value is 1°. The sample was exposed to 
an ion fluence of 2 x 10^^ ions ni~^, with an average ion 
flux of 3 X 10^^ ions ni~^ s~^. The flux is specified here 
for the sample plane. It was controlled by a Faraday cup 
movable to the sample position. After ion exposure, the 
samples were imaged in situ by STM. 

A representative image is shown in Fig. [TJ Superficial 
naked-eye inspection already suggests the existence of a 
long-wave corrugation which does not seem particularly 
ordered but does have a preferred orientation, that coin- 
cides with the projection of the ion-beam onto the sub- 
strate plane. Actually, this is taken into account in the 
analysis of the topography. Specifically, the ion beam 
comes from the lower left corner with an angle of 35° to 
the X axis. After rotating the images so that the ion- 
beam projection coincides with the x (horizontal) axis, 
rectangular region have been cut out from the original 
figures and provide the data for our analysis, see jis} . 

Projections of the two dimensional PSD along the 
and ky axes are presented in Fig. [21 These plots have been 
obtained after averaging the S'(k) functions obtained for 
different window sizes, as detailed in the Supplemental 
Material, to improve statisticsJ^ Note that power-law 
behavior takes place for length-scales that are larger than 
approximately 6 nm, while for shorter distances (larger 
corresponding values of k) different behavior is obtained. 
Thus, we will restrict our analysis in this case to the 
small k range, thus indeed facing the properties of the 
long- wavelength corrugation just mentioned. 

By performing a fit over the small wave-vector region 
in Fig. [21 we estimate the values of the momentum-space 
roughness exponents to be 

2(5^ 2.66 ± 0.02, 25„ = 1.80 ± 0.02, 

(16) 

C = 1.48 ±0.02, 

where error bars come from statistical uncertainty in the 
fits. Hence, strong anisotropy occurs, C 1- The corre- 
sponding exponents characterizing the power-law decays 




k (nm' ) 

FIG. 2. One-dimensional projections S{k, 0) (green circles, 
left axis) and S{0, k) (blue squares, right axis) of the two- 
dimensional PSD for the contamination-free sample, averaged 
over 8 windows. The solid red and dashed black lines are 
fits obtained by least squares, the resulting exponents being 
2ax — 2.66 and 2ay ~ 1.80. These lines correspond to the 
exact stationary behavior of Eq. (|19p for = 178, Vy — 
0.99, and D = 1 with nm and s for space and time units, 
respectively. 

of the ID PSD functions should be 

2a^ + 1 = 1.18 ±0.02, 2a3, + 1 = 1.12 ±0.02. (17) 

Using these values in Eqs. © and pH)) . the scaling be- 
havior of Sy{ky) extracted from our Ansatz agrees well 
with experimental data as shown in Fig. [31 while such 
an agreement is not reached in the case of S^lk^). In 
order to understand this disagreement, note first that in 
the present system both real-space exponents are quite 
close to being effectively zero, in view of the statisti- 
cal fiuctuations in the data. In the context of criti- 
cal phenomena and power-law behavior this would be 
associated with logarithmic (rather than proper power- 
law) behavior for the real space correlation functions 
Incidentally, when issues like this arise, the assessment 
of kinetic roughening behavior is usually more clear-cut 
for Fourier-space observables than for real-space correla- 
tion functions, see [l^ and references therein. On the 
other hand, as we can see in Fig. [31 the scaling behavior 

predicted by the simple power law Sx{kx) ~ /c^; 
has not been reached by the experimental data, espe- 
cially for the smallest wave- vector values. An analogous 
inconsistency between the scaling behavior of the two- 
dimensional PSD and that of the ID PSDs has also been 
observed in continuum models of SA surfaces.— The key 
to understanding the discrepancy is the fact that, for the 
small numerical values obtained for ax,y, integrals such as 
([6]) converge slowly to their asymptotic properties, in the 
sense that non-negligible corrections occur to the lead- 
ing terms ^"(^"^.I'+i) in a small kx,y expansion. In other 
words, one needs to rather compute the full integrals ^ 
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k (nm' ) 

FIG. 3. PSD of one-dimensional cuts Sx{k) (green circles, 
left axis) and Sy{k) (blue squares, right axis) of the 6 = 81° 
sample, averaged over 8 windows. The solid magenta line and 
the dot-dashed brown line correspond to Eqs. (|9]) and (|10|l . 
respectively, with values given by Eq. (|17|l . and are not fits 
to the data. Rather, they have been computed using p8p for 
Eq. gH) with = 178, I'y = 0.99, and D ^ 1, with nm and 
s for space and time units, respectively. 

that, for a finite space resolution, read 

|.7r/a 

Sx,y{kx,y) I dky X S{kx , ky) , (18) 
Jo 

where a is the smahest length-scale in the region ex- 
plored. Even in such a case, convergence to the asymp- 
totic behavior may occur only at length scales that may 
not be easily accessible in an experiment. 

In the present case, based on results from we can 
write down a continuum model with the same scaling ex- 
ponents as measured in , for which S A behavior holds 
rigorously, and analyze for it the behavior of the ID PSD 
functions. Thus, as a proxy for the scaling behavior of 
the present experimental data, we can write down the 
following linear evolution equation for the Fourier com- 
ponents of the height field, ft.(k, i), 

dth = -(i/,|fc.p-66 + z.JfcJi-80)/i + f,. (19) 

Note, this equation has a similar structure to the one Eq. 
([5]) would have in Fourier space, only with values of the 
exponents that are those from p^ . instead of 2ax — 2 
and 2cey = 4, as would be the case for the latter. Eq. 
(1191) has a complicated expression in real space, and can 
be interpreted as a model of anisotropic non-local surface 
relaxation.^- It can be readily shown^^ that, in the long 
time limit, Eq. ([T^ follows the scaling Ansatz ^ exactly 
with exponents given by (jl6p . Actually, the power-law 
fits in Fig.[2]are the exact stationary behavior of Eq. ([T9|. 
provided we use Vx = 178, Vy = 0.99, and D — 1, with 
nm and s for space and time units, respectively. Once this 
is the case, the behavior of the ID PSD functions for Eq. 
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FIG. 4. AFM top view of the morphology of Si(lOO) irradiated 
for 8 h with 40 keV Ar''' ions and concurrent Fe codeposition. 
Image size is 6 x 6 fim^. 

([T9l) follows without the possibility to further fits, and is 
shown in Fig. |3l As we see, there is good quantitative 
agreement with experimental data for Sy{ky), while in 
the case of Sx{kx) the experimental data overshoot the 
theoretical curve at small kx values. For a full class of 
models related with Eq. (IT5|) this is an indication that 
the true asymptotic behavior is reached at still larger 
distanceSfl^ and we beUeve this is the case for the present 
set of data. Note, we are not claiming that Eq. (fT9|) pro- 
vides the description of the full dynamics of the latter. 
We do believe it shares with it the same asymptotic be- 
havior, and this is possible even if the "true" dynamical 
equation for this system is a different, even nonlinear, 
onej^ An interesting physical conclusion for the model- 
ing of IBS problems stems from the fact that Eq. 
preserves the average value of the height field along the 
system evolution, since {h{0,t)) = {h{t)) = const. This 
suggests in particular that erosive mechanisms^^ that 
do not preserve this quantity, are (numerically) negligible 
under the present experimental conditions, which agrees 
with recent experimental results for similar systems, see 
G-g- [ill references therein. 



B. Morphologically unstable case 

The irradiation experiments in the morphologically un- 
stable case were performed with a 40 keV Ar+ beam ex- 
tracted from a Danfysik 1090 ion implanter with a base 
pressure of 5 x 10^ mbar, similar to ^l[. The ions im- 
pinged on the single-crystal Si(lOO) targets (1x1 cm^) 
at 60° ± 5° with respect to the surface normal with a 
current density of 18 /iA cm~^ in the sample plane. A 
steel plate (1.5 mm high) placed adjacent to the Si target 
acted simultaneously as Fe source and sample holder. To 
obtain homogeneous irradiation, the focused beam was 
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k (nm' ) 

FIG. 5. One-dimensional projections S{k, 0) (red circles, 
left axis) and S{0, k) (blue squares, right axis) of the two- 
dimensional PSD of the samples illustrated in Fig. 3] after an 
average over 7 windows. The black dashed line and the solid 
magenta line are fits obtained by least squares, the resulting 
exponents being 2ax — 5.81 and 2ay ~ 3.74. 




k (nm' ) 

FIG. 6. PSD of one-dimensional cuts Sx{k) (red circles, left 
axis) and Sy{k) (blue squares, right axis) of the samples il- 
lustrated in Fig. 31 after an average over 7 windows. The 
black dashed line and the solid magenta line correspond to 
Eqs. ((9]) and (|10p . respectively, with exponents values given 
by Eq. (gTJ. 



scanned with a magnetic x — y sweeping system in such 
way that both the Si surface and the steel target were 
bombarded. The irradiation time was 8 hours. The re- 
sulting surface morphology was imaged ex- situ by AFM 
operating in the dynamic mode with a Nanoscope Ilia 
equipment (Veeco®). Silicon cantilevers, with a nominal 
radius r of 8nm and opening angle smaller than 52°, were 
employed. A representative topography is shown in Fig. 
m The projection of the ion beam coincides with the x 
(horizontal) axis in the figure. 

Although from the physical point of view the present 
system is perhaps more complex (due to the not well 
understood role of impurities in the nanopatterning pro- 
cess), the occurrence of SA seems in principle clear-cut. 
We employ the same procedure as above. After comput- 
ing the two-dimensional PSD, we show its two projections 
S'(fca;,0) and S{Q,ky) in Fig. [5] Both plots appearing 
in the latter figure immediately suggest the existence of 
a characteristic length-scale around 800-900 nm, above 
which no proper scaling takes place. We associate this 
scale with the wavelength of the pattern that develops 
under these experimental conditions. Nevertheless, at 
smaller scales (larger k values) power-law behavior ap- 
pears that is different in the and ky directions, as 
described by exponent values 

2(5^ = 5.81 ±0.06, 2aj, = 3.74 ± 0.2, 

(20) 

C = 1.55 ±0.01, 

that are obtained by least squares. Note, C 1- Thus, 
according to Eqs. (IT^ and p^ . exponents characterizing 
the power-law decay of the ID PSD functions should be 

2q;2; + 1 = 4.26 ±0.14, 2ay + l = 3.10±0.2. (21) 



Using these values in Eqs. ^ and (ITUl) . the agree- 
ment with the experimental data obtained for the one- 
dimensional PSDs for the corresponding range in k is 
very good, as shown in Fig. [HI 

The appearance of scaling behavior at length scales 
below the one associated with the morphological insta- 
bility suggests as a cause the influence of the average 
shape associated with it, rather than kinetic roughen- 
ing, similar to the case of mound formation in epitaxial 
growth on high symmetry surfaceSii In order to explore 
this possibility, we have generated artificial topographies 
made up by an array of semi-ellipsoids that are placed 
off their ordered positions by a random error, or a dis- 
ordered array of faceted featuresJ^ In both cases, while 
the interplay between disorder and the smooth shape of 
the basic motif (semi-cllipsoid or sawtooth) does give rise 
to power law behavior of 2D and ID PSD functions be- 
low the corresponding characteristic sizes, the scaling ex- 
ponents do not fulfill the relations expected within our 
scaling Ansatz. Rather, scaling in these cases seems re- 
lated with the smooth geometry of the basic patternj^^ 
combined with fluctuations in its space arrangement. 

We thus suggest strongly anisotropic kinetic roughen- 
ing at submicron scales (above which scale invariance is 
lost) to account for the present experimental system fea- 
tures. To some extent, such a behavior is complemen- 
tary to the one that is typical of e.g. the KS system,— in 
which a small-scale pattern becomes disordered at suffi- 
ciently large scales at which kinetic roughening occurs 
Although in terms of the latter class of systems it may 
look somewhat peculiar, the behavior found for the ex- 
periments just analyzed is again readily obtained in ap- 
propriate simple model systems. Thus, one can general- 
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ize e.g. Eq. ([5]) to 

dth = ~fih + Vxdlh - Vydyh + ■q{x, y, t), (22) 

where the Unear term with coefficient fj, > accounts 
for a physical mechanism (e.g. a wetting layer) that fa- 
vors a specific value of the height (say, h = 0). For 
Eq. (1221) . at i — > cxD one can exactly obtain S{kx,ky) ~ 
{fj, + Vxk'^ + i^yky)~^ . This implies that S{kx,ky) does not 
scale with fc^;,!, (but, rather, becomes a constant) for the 
largest length-scales (smallest k^^y values), thus break- 
ing scale invariance at such large scales. However, ki- 
netic roughening precisely as in Eq. does still occur 
for (large) k^^y values whose contribution to S{kx,ky) 
dominates over fi. While we are not saying that Eq. ((22|) 
describes the experiments just analyzed, we believe the 
latter may correspond to a situation that is qualitatively 
captured by this example. Currently there is a large ef- 
fort in pursuing a more quantitative description of IBS 
experiments in the presence of metallic contamination 
that takes into account these type of effects, see [1^ and 
references therein. 



III. DISCUSSION AND CONCLUSIONS 

We have verified a scaling Ansatz for strongly 
anisotropic kinetic roughening properties for two experi- 
mental systems. First of all, our results provide a consis- 
tent experimental assessment of such type of behavior, 
which calls for the development of theoretical models 
that can describe in detail the experimental properties 
observed. This would add to an improved knowledge on 
the phenomena driving the specific systems considered, 
as well as to our understanding of driven systems in gen- 
eral. From a more practical point of view, our work il- 
lustrates several forms of occurrence of SA, both in the 
absence and presence of morphological instabilities, and 
at (sufficiently large) scales that remain either above or 
below a privileged (pattern) size. 

Generally speaking, even in absence of typical length- 
scales, it is important to keep in mind that our scaling 
Ansatz applies only in a hydrodynamic limit. There- 
fore, in actual physical systems, effects due to finite-size 
and/or finite space resolution of measurement techniques 
could hinder a clear-cut scaling of the correlation func- 
tions, thus preventing the observation of a consistent be- 
havior between the one and two-dimensional PSDs. As 



we have seen, this "slow convergence" issue appears to be 
more pronounced for integrated quantities, such as the 
one-dimensional PSDs [see formula ([SJ], and has been 
studied in more detail in Tt'] by means of analytical 
arguments as well as numerical simulations. Thus, the 
two-dimensional PSD, conveniently averaged over several 
samples (or windows of the same sample), appears to be 
a more reliable observable in order to asses anisotropic 
scaling. However, in experiments it is seldom possible to 
have a sufficient number of samples over which to per- 
form such an average, and the two-dimensional PSD fre- 
quently appears too noisy to extract reliable values of 
the exponents. The one-dimensional PSDs have a higher 
signal-to- noise ratio due to the fact that, for each win- 
dow, they are averaged over all the rows or columns of 
the matrix. For this reason, despite the fact that they are 
more sensitive to finite size effects, such quantities do in- 
deed provide relevant information about surface features. 

We can hypothesize that this type of complications 
may have hindered a more frequent observation of 
anisotropic scaling in the literature. We hope that, once 
the analysis has been clarified, the identification of this 
challenging type of behavior becomes simplified and we 
can thus understand it better. Perhaps an analogy can 
be drawn with the case of anomalous scaling. This is 
a type of kinetic roughening behavior that was usually 
associated with large values of the roughness exponent, 
that induced many wrong assessments of scaling behav- 
ior both in experiments and models.— Incidentally, such 
"anomalous" values usually introduced slow convergence 
properties in correlation functions. Once anomalous scal- 
ing was identified and systematized^"— it has been in- 
deed assessed in different thin-film experimental systems, 
from e.g. electrochemical depositions^ to chemical vapor 
deposition,— metal dissolution^^ etc. We can only hope 
for a parallel development regarding strongly anisotropic 
surface kinetic roughening in the near future. 
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